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- fonctionnement d’un turboreacteur 

- organes d’un turboreacteur 



Le dimensionnement 




Rappel : principe de fonctionnement d’un turboreacteur 
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Rappel : principe de fonctionnement d’un turboreacteur 



Compresseur Turbine 

Soufflante haute pression haute pression 




Axe de turbine 
basse pression 



Tuyere 



Axe de turbine 1 
haute pression 



Compresseur 
basse pression 



Chambre de 
combustion 



Turbine 
basse pression 




Organes structuraux principaux 






Organes structuraux principaux 



Rotor BP 



Rotor HP 






Stator 




Organes structuraux principaux : rotor 




1 8 3 4 5 6 




(a) Rotor : 1 — paLier avant du compresseur 
(rouleaux), 2 — rotor du compresseur. 3 — pa- 
lier arriere du compresseur (billes), 4 - arbre de 
la turbine, 5 — palier de la turbine, 6 - rotor de 
la turbine 




TWIN -SPOOL COMPRESSOR 



MOUNTING FLANGE 



H R COMPRESSOR 



H.R SHAFT 

DRIVE FROM TURBINE 








Organes structuraux principaux : stator 



1 2 3 S 6 7 




(b) Stator : 1 — carter d ’entree du compresseur. 
2 — carter du compresseur. 3 — carter axriere 
du compresseur, 4.5 — carters de la chambre de 
combustion, 6 - carter d’ entree de la turbine, 7 
— carter de la turbine 







Preliminaires 

- fonctionnement d’un turboreacteur 

- organesd'un turboreacteur 



Le dimensionnement 




Introduction 



Un constat 



• si la performance est surtout aero... 




c ar act£ r istiqu es 




■ Pau&see mait, an d^collage (kN) 


42$ 


* T auH de dilution 




■ k poL&s^e nominate 
mairrtcnue fQ 


30 


* Pou5s4e max. en montee (k[\|) 
35 000 ft - Mach 0,84 ISA 
Moteur avion rse 




* T&ux d e dcimpiesiion general 
pour la monte-e max. 


40 


- Longueur (mm) 


7 289 


■ Diametre soufflante (mm) 


3 403 


* Applications 


B777-20Q 

B/77-200ER 



• la defaillance est souvent structurale... 




Rupture par HCF sur banc d’essai de turbine a gaz 




Rupture par HCF sur LPT3, Boeing 747-338, 2004 



www.snecma.com / M. Jocker, these KTH, 2002 - Rolls Royce pic / www.atsb.gov.au 










Introduction 



Historique 



Quelques chiffres 



“I am going to make a gas turbine before I quit and it is going 
to have higher efficiency than the reciprocating gas engine, 
and it is going to run without vibration ...” (Glenn Warren, 
1914). 

‘‘On axial compressor blade vibration we have gone very 
cautiously and I believe have not met the problem in all its 
bleakness; the dark days are still ahead” (Hayne Constant, 
1947). 

‘‘Suppression of vibratory fatigue is the design challenge” 
(Clarence Danforth, 1967). 

“Failure of fan blades on the B-l’s General Electric engines 
is the suspected cause in both incidents . . .” ( Los Angeles 
Times, 1991). 

“. . . [T]he engine failures are a continuation of fourth stage 
turbine blade cracking problems that first surfaced in the Pratt 
and Whitney F100-229 powerplant ...” (Aviation Week & 
Space Technology, 1994). 

“Engine high cycle fatigue is the Air Force’s #1 readiness 
issue” (John Deutch, 1995). 



Modes de defaillance des 
turboreacteurs d’avion 



FOOflKO,** 




- 90% of the potential HCF problems are covered during development testing the 10% remaining problems account 
for nearly 30% of the total development cost and cause over 25% of all engine distress events (El-Aini et Al, 1997) 

-every new development programme for jet engines has about 2.5 serious HCF problem (Kielb, 1998) 

- US Air Force and Navy estimate that cost to them approximately US$400 Million per year (Garrison, 2001) 

A. V. Srinivasan, JEGTP, 1997 / Fransson, FUTURE project, 2007 / Epstein, Aeromechanics lectures, 2004 






Introduction 



Contexte industriel 



Evolution des turbomachines modernes : 

- diminution de la masse 

- augmentation de la charge 

- reduction des dimensions axiales 



Choix technologiques : pieces monobloc 

- desaccordage inter-aubes 

- amortissement faible 




Consequences : 

- reduction des marges au flottement 

- niveaux dynamiques plus eleves et plus 
difficiles a predire 

- augmentation des interactions entre les roues 




Figure 1.1-2: Trends in turbine blade design (right) and example of a modern slender (bottom left) 

and old (top left) turbine rotor. 



FUTURE project (http://www.future-project.eu), pictures from Volvo Aero (http://www.volvoaero.com) 






Introduction 



Phenomenes et modeles pour le dimensionnement d'aubages 




regime 

[runsitoire 




nature 

chfLutique 





decorr£li£s avec la 
Vitesse 4c rotation 



■ variation dc la 
vitcssc dc 
rotation 



- couches Limiles 

- Ullages 
-jet& h cie , 



- touibi lions dc 
Von Kantian 



- pompage 

- fluTtcmcnl 



• dans la suite : focus sur 
I’estimation de la DDV en 
fonctionnement « normal » 

mais de nombreux outils sont 
aussi utiles aux limites des 
domaines d’exploitation 



Natural Frequencies 
Mode shapes 





S. Callot, these ECL, 2002 



J. S. Rao, INCAST, 2008 



4 





Introduction 



Dimensionner des aubages = un processus pluri-disciplinaire 




M. Jocker, these KTH, 2002 






Plan 



A. Pre-design aero 

B. Analyse du comportement sous chargement statique 

C. Estimation du risque vibratoire 

D. Tests sur prototypes 





A. Pre-design aero 



Exemple de la conception d’une RM destinee a un etage de compresseur 

Cahier des charges 



• environnement « moteur » / banc 



Maximum power 


800 kW 


Rotational direction viewed FLA 


Clockwise 


Maximum torque 


340 Mm 


Maximum shaft speed 


20300 rpm 





^performances ciblees 



Design speed 


1 8000 rpm 


Stall margin 


> 20% 


Flow coefficient 


-* 0.6 


Reduced frequency 


~ 0.25 


Target frequency at ADR for, 1 F 


320 Hz 


Flutter limit 


Hadp < n HuL1t= , < n sla „ 


Available max torque sufficient to reach 


surge 

105% speed 



* contraintes mecaniques 



Tip clearance 


0.5 mm 


Maximum allowable stress 


>60% of yield stress 


Disc burst margin 


1.4 



FUTURE project (http://www.future-project.eu), from Volvo Aero (http://www.volvoaero.com) 
and PCA Engineers Limited (http://www.pcaeng.co.uk), rigs’s specifications in the first chart 
from TUD - Technischen Universitat Darmstadt (http://www.tu-darmstadt.de/) 
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A. Pre-design aero 



Exemple de la conception d’une RM destinee a un etage de compresseur 



Proposition et evaluation d’un profil d’ALJBE 

• definition d’un profil 





le redresseur est suppose connu 
nombre d’aubes du blisk choisi 



etude aero stationnaire 



Total, pressure 



Entropy - 'JrjM, span 




I 




Figure 6: Steady NS computation results 



prediction des performances 
de I’etage 





[kg/s] 


Max [], r H 


SM [%] 


iJr.pv* [%] 


’h.ADF [%1 


13.9 


1.476 


1*3 i 


91.3 


90.8 



Wilson and Korakianitis, Cascade design parameters, 1 998 / J.-P. Lame / Seinturier etal, GT 2004 / Leboeuf etal, ICAST 201 0 





Plan 



A. Pre-design aero 

B. Analyse du comportement sous chargement statique 

C. Estimation du risque vibratoire 

D. Tests sur prototypes 




B. Analyse du comportement sous chargement statique 



Exemple de la conception d’une RM destinee a un etage de compresseur 



Proposition et evaluation d’un profil de DISQUE 



STRUCTURAL UNSTEADY AERODYNAMICS STEAQY AgRQPYHAMilCS 




• I’analyse aero stationnaire fournit les chargements pour calculer 

• les niveaux de contrainte moyens subis par I’aube et le disque 

• les deplacements statiques en fonctionnement 

1/ jeu au carter : compromis pertes aero / risque de contact (utilisation de revetements abradables) 
21 modification des profils sous charge : reverse design pour obtenir le profil a fabriquer ! 

• processus iteratif / optimisation 



• contraintes mecaniques 



Tip clearance 


0.5 mm 


Maximum allowable stress 


>60% of yield stress 


Disc burst margin 


1.4 



Schematic of forced response prediction system [Hilbert et al. 1997] Projet FUTURE, resultats intermediates 






B. Analyse du comportement sous chargement statique 



Exemple de la conception d’une RM destinee a un etage de compresseur 



Proposition et evaluation d’un profil de DISQUE 




FUTURE project (http ://www.futu re- 
project. eu), analysis performed by 
PCA Engineers Limited 
(http ://www. pcaeng.co.uk) 



CftfntMulim* EipaKs(*H Eihmt 

3000 

2000^ 

A. Sternschuss, these 
I mm I. ECP, 2010 

E 

a 

• ambiance thermique dans un moteur 




CFM 56-5A 


Compresseur BP 


Compresseur HP 


Turbine HP 


Turbine BP 


Aubes stator 


Acier 


Acier -> nickel 


Nickel 


Nickel 


Aubes rotor 


Titane 


Titane -> nickel 


Nickel 


Nickel 



• exemple d’utilisations de materiaux 






B. Analyse du comportement sous chargement statique 

Exemple de la conception d’une RM destinee a un etage de compresseur 

Proposition et evaluation d’un profil de DISQUE 

* effets centrifuges : une ailette de 0.2kg situee a un rayon de 0.3m de I’axe moteur 
tournant a 10000 rpm subit une force centrifuge equivalente a 9t 



10000 



1000 7 



2. Resistance-Densite 



' Ceramiques^'3; :i 7,>, 
technique 

•J5i X7 /**, 





Aube RTM tissee 3D 



materiau retenu : alliage de titane Ti6-4 



Young’s modulus- 


110-114 GPa 


Poisson's ratio 


0.31 


Density 


4420 kg/m J 


Ultimate tensile stress 


395-930 MPa 


Yield stress 


825 -859 MPa 


Ultimate strain 


6-10% 



Densite p(Mg/m 3 ) 



^j| 



Ashby, Choix des materiaux, Dunod / www.snecma.com 



Projet FUTURE, resultats intermediaires 




B. Analyse du comportement sous chargement statique 

Exemple de la conception d’une RM destinee a un etage de compresseur 

Proposition et evaluation d’un profil de DISQUE, verification de la tenue statique de I’aube 

• parts relatives du centrifuge et de I’aero moyen Encyclopedia of vibrations, 
sur le niveau de contrainte en pied d’aube Elsevier, 2002 



Steam turbine 




Biaried disk: last stage, low pressure 










Power 

(MW) 


Speed 

(rpm) 


Power 

(MW) 


No. of 
biad&s 


Blade length 
(mm) 


Blade mass 

{kg) 


Base centrifugal 
force (W) 


Maximum static base strain 
(Nm -2 ) 
















Centrifugal force 


Torque 


360 (fossil) 


3600 


17 


78 


865 


23 


2 x Iff 5 


i.5 x icy 3 


0.12 x 10 3 


1000 (nuclear) 1500 


19.6 


136 


1220 


39 


1.22 x 10® 


0.276 x 10 a 


0.123 x 10° 



reponse statique au centrifuge et a I’aero moyen 



critere en contrainte (VM) 



(calcul en NL geometrique) 



odg : 350 MPa 




23:3B PP Standard T*nr ETO5 



hhj-1, CF at ISBCbpm 
"t 5r r.rrt* - |(Ul 
•Jot: 

hz vj t: U rx-fj-n-i jT'iin Scab- Farter 



+ video blisk 



critere en deplacements (vue radiale) 




- effet de devrillage (odg : 1mm) 
radial + correction thermique (odg : 0.5 mm) 



FUTURE project (http://www.future-project.eu), analysis performed by PCA Engineers Limited (http://www.pcaeng.co.uk) 






Plan 



A. Pre-design aero 

B. Analyse du comportement sous chargement statique 

C. Estimation du risque vibratoire 

D. Tests sur prototypes 




C. Estimation du risque vibratoire / Demarche 





Techniques de I’lngenieur, BM 4 175 




Plan 



A. Pre-design aero 

B. Analyse du comportement sous chargement statique 

C. Estimation du risque vibratoire 

I. Calcul des modes propres 

II. Determination des frequences d’excitation 

III. Recherche des coincidences 

IV. Calcul des contraintes dynamiques 

V. Trace du diagramme de Haigh 

D. Tests sur prototypes 




C.IV. Estimation du risque vibratoire / Calcul des modes propres 



Demarches 



• Hypotheses 

- calculs sous vide d’air 

- etude d’une roue isolee 

• Analyse modale d’aubes isolees 

- prise en compte des effets de la rotation 

• Analyse modale en symetrie cyclique 

- pour des assemblages accordes 

- prise en compte des « comportements de groupe » et des couplages aubes-disque 

- etude d’un seul secteur : aube + disque / conditions aux limites imposees periodiques 

•Analyse modale complete 

- pour des assemblages desaccordes 

- techniques de sous-structuration et de reduction (Craig-Bampton...) 

- tres couteux 





C.IV. Estimation du risque vibratoire / Calcul des modes propres 



Aube isolee 



Effets centrifuges / gyroscopiques / spin-softening (cf cours J.-P. LaTne) 
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+ video LDV 



FUTURE project (http://www.future-project.eu), analysis performed by PCA Engineers Limited (http://www.pcaeng.co.uk) 






C.IV. Estimation du risque vibratoire / Calcul des modes propres 



Etude en symetrie cyclique 

Effets de groupe / couplage aubes-disque (cf cours F. Thouverez) 



Le nombre de diametres nodaux est egal au nombre de maxima de la deformee sur la circonference 

Zone ou seul I’effet de groupe 




FUTURE project (http://www.future-project.eu), analysis performed by Ecole Centrale de Lyon (http://www.ec-lyon.fr) 






C.IV. Estimation du risque vibratoire / Calcul des modes propres 



Etude en symetrie cyclique 



Representation 3D des modes propres 





C.IV. Estimation du risque vibratoire / Calcul des modes propres 



Etude globale de la roue aubagee 

Necessity de modeliser le mistuning 




Engfna Order 

Enveloppe des amplifications de la reponse forcee 
imputables au mistuning en fonction de I’ordre des 
excitations (EO = nombre de cycles d’excitation par tour) 



• simulation 

• modeles probabilistes parametriques 

• identification des motifs de mistuning 

• pour le recalage 




Figure 3: Mistuning pattern 



blade 



* interet du mistuning intentionnel 



Castanier & Pierre, JPP 2006 - F. Pichot, these ECL, 2002 






C.IV. Estimation du risque vibratoire / Calcul des modes propres 




Comment tenir compte, en simulation, du couplage inter-etages ? 




Figure 1.12: History of multi-stage computations at Snecma: (a) Fourier axisymmetric* (b) 3D 
slice, (c) 3D cyclic symmetry, (d) 360° shell/volume model and (e) multi-stage cyclic symmetry 



Couplage des bases cycliques 



A. Sternschuss, these ECP, 2010 



Laxalde, these ECL, 2007 







Plan 
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I. Calcul des modes propres 
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C.IV. Estimation du risque vibratoire / Frequences d’excitation 



Demarche 



• Notion d’Engine Order (EO) : 
nombre de perturbations du flux / tour 

- vu de I’organe considere 

- vision harmonique en 6 

- excitations synchronises avec la vitesse de rotation 

• LEO (1 a 8 harmoniques) 

FREQUENCY 

- distorsion du flux d’entree 

- bras de carter 

- prelevements d’air 

- injecteurs (turbine) 

- balourd (roues fixes) 



• HEO (depend du nombre d’aubes des roues voisines) 




- RM : sillages de 2 stators amont et 1 aval 

- RD : sillages de 2 RM amont et 1 aval 

- attention : la « transformee de Fourier du sillage » en 6 fait 
que les multiples entiers des nombres de sillages doivent 
aussi etre pris en compte (2 e harmonique surtout) 



LNG1I 




Jay etal, AGARD, 1 988 / Leboeuf etal ICAS 2010 / R Payer, 201 0 
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C.IV. Estimation du risque vibratoire / Coincidences 



Demarche 



• Conditions necessaires pour la resonance d’une roue aubagee 

• une coincidence frequentielle : 

frequence de I’excitation = frequence du mode 
-> diagramme de Campbell 

• ET une coincidence spatiale : 

nombre d’obstacles de I’excitation = nombre de diametres du mode propre 
(modulo le nombre d’aubes) 

-> diagramme zig-zag 





C.IV. Estimation du risque vibratoire / Coincidences 



Demarche 

Conditions de coincidence spatiale 

• le diagramme zig-zag 

• permet de determiner le nombre de diametres nodaux caracteristique des modes excitables par un 
nombre d’obstacles donne (sillages ou premiers harmoniques) 

• si la structure excitee etait axisymetrique (disque plein), une excitation tournante periodique 
exciterait les modes dont le nombre de diametres nodaux est egal au nombre d’excitations par tour 

• pour une roue aubagee ce sont les aubes qui sont excitees ; d’ou un nombre de diametres apparent 
borne par N/2 



Engine order 

(nombre d’obstacles / tour de I’excitation) 




Nombre de 
diametres du 
mode excite 



Adapte d’un document de formation Snecma 






C.IV. Estimation du risque vibratoire / Coincidences 



Demarche 

Condition de coincidence spatiale 

• illustration 

• des fonctions de forgage a differents EO 
excitent un meme mode a 2 diametres 

• N = 20 aubes sur le disque 



N=20 D=2 Dex=alpha*N-D or alpha*N+D 




Engine order 

(nombre d’obstacles / tour de 
I’excitation) 




(N-1J/2 





Nombre de 
diametres du 
mode excite 



• principe de la demonstration 
dans le cas general 

• trouver les conditions d’extremalite des 
forces general isees associant le mode et 
la forme d’excitation consideres 



F. Moyroud, these KTH 1998 






C.IV. Estimation du risque vibratoire / Coincidences 



Demarche 



Representation 3D du forgage et des conditions de coincidence 
spatiale (pour le premier segment du zig-zag, EO=ND) 






C.IV. Estimation du risque vibratoire / Coincidences 



Demarche 



Exploitation du diagramme zig-zag (sur son premier 
segment EO<ND/2), a vitesse de rotation donnee, EO=20 



7000 



Frequence = 20xCO 




• il y a probleme 

• plusieurs solutions sont possibles : 

• changer le nb d’aubes sur le redresseur amont 

• modifier la raideur des aubes de la roue mobile 

• modifier la raideur du disque 

• changer la vitesse de fonctionnement nominal 



20 ND excitation line 



0 10 20 30 40 50 

NODAL DIAMETERS (ND) 




wm 



SAFE diagram, M. P. Singh etal, Dresser-Rand Co website 
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C.IV. Estimation du risque vibratoire / Calcul des contraintes dyn. 



Demarche 



• Beaucoup d’inconnues en debut de conception 

• complexity du comportement de la structure (mistuning, non-linearites) 

• estimation de I’amortissement total: combinaison des amortissements structural et aeroelastique 

• un calcul precis de la fonction de forgage aerodynamique est requis 

• Points-clefs dans le calcul de la reponse forcee 

• quantification des dissipations structurales 

• I’ aerodynamique instationnaire / ingredients du modele 

• le couplage fluide-structure 



Schematic of forced response prediction system [Hilbert et al. 1997] 




STRUCTURAL UNSTEADY AERODYNAMICS STEADY AERODYNAMICS 






C.IV. Estimation du risque vibratoire / Calcul des contraintes dyn. 



Estimation des niveaux cTamortissement 

Importance de I'amortissement en dynamique 



Amplification dynamique de la FRF 
d’un oscillateur a 1 ddl avec un 
taux d’amortissement % 

0(11) ~ , — ; f2 = 

V(i -n“) 2 + 4^n* 




• I’amplitude de la reponse a la resonance est controlee par I’amortissement 

• pour £«1 I’amplification est inversement proportionnelle a 2, 

• un soin particulier doit etre apporte a sa modelisation 

• importance d’un recalage par mesures experimentales 



C. Soize. Dynamique des structures. Ellipses 





C.IV. Estimation du risque vibratoire / Calcul des contraintes dyn. 



Estimation des niveaux d’amortissement 

Fluide vs structure 



FREQUENCE en hertz 











TORS 


Q*i 






























































, f .. 
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r * k 







aoL 




sans VENT - Amortisae-nrienl macefiiquv 




ay PC VENT V = 131 mfa 
Amitr+isscment aarodynamiqua 





• amortissement structural du au materiau : ^<0.001 si monobloc (hysteretique ou visqueux) 

- EXPERIMENTAL 

• amortissement structuraldu a d’autres sources de dissipation (frottement) 

- EXPERIMENTAL + CALCUL 

• amortissement aero predominant devant amortissement materiau (visqueux, parfois « negatif » si instability) 

- CALCUL + EXPERIMENTAL 

• necessity de dispositifs ad hoc, en particulier pour le monobloc 



M. Pluviose. Ingenierie des turbomachines. Ellipses 2003 






C.IV. Estimation du risque vibratoire / Calcul des contraintes dyn. 



Estimation des niveaux cTamortissement 



Solutions technologiques pour augmenter I'amortissement 



• Dispositifs lineaires 

• revetements absorbants (visco-elastiques, elastomeres, silicones) sur aubes ou aux interfaces 

• fluides dans des aubes creuses (experimental) 

• Dispositifs non-lineaires 

• par frottement : en pied d’aube, par nageoires, talons, cables, frotteurs (annulaires = joncs, plaques, doigts...) 

• par chocs (experimental) 

• par controle passif ou semi-actif (experimental) 

• Autres : mistuning intentionnel 




S. Nacivet, these ECL/ D. Laxalde, these ECL 






C.IV. Estimation du risque vibratoire / Calcul des contraintes dyn. 



Estimation des niveaux cTamortissement 



Solutions technologiques pour augmenter I'amortissement 

• Dispositifs lineaires 

• revetements absorbants (visco-elastiques, elastomeres, silicones) sur aubes ou aux interfaces 

• fluides dans des aubes creuses (experimental) 

• Dispositifs non-lineaires 

• par frottement : en pied d’aube, par nageoires, talons, cables, frotteurs (annulaires = joncs, plaques, doigts...) 

• par chocs (experimental) 

• par controle passif ou semi-actif (experimental) 

• Autres : mistuning intentionnel 




Fig . 1 . 17 : Amortisseur par frottement de type « doigts » pour les disques de turbines monoblocs ; Brevet 
Boeing US 6.375,428 B1 



D. Laxalde, these ECL / Rao, Turbomachine blade vibration, 1 991 





C.IV. Estimation du risque vibratoire / Calcul des contraintes dyn. 



Estimation des niveaux cTamortissement 

Solutions technologiques pour augmenter I'amortissement 

• Dispositifs lineaires 

• revetements absorbants (visco-elastiques, elastomeres, silicones) sur aubes ou aux interfaces 

• fluides dans des aubes creuses (experimental) 

• Dispositifs non-lineaires 

• par frottement : en pied d’aube, par nageoires, talons, cables, frotteurs (annulaires = joncs, plaques, doigts...) 

• par chocs (experimental) 

• par controle passif ou semi-actif (experimental) 

• Autres : mistuning intentionnel 
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Estimation des niveaux d’amortissement 

Calcul des reponses non-lineaires: exemple du frottement 




a I’interieur des solides 



equilibre interne 
comportement (elastique 

lineaire) 



aux interfaces 



non penetration (usure = 
jeu variable) frottement 
tangentiel (Coulomb) 
loi d’usure (Achard) 



sur les autres frontieres 



deplacements imposes 
excitation periodique 



M U (l) + CU {,) + K U {,) + F c (,) (U {,) .U (, \lV {,) ) = f£> 

i 

• expression de I’equation du mouvement dans le domaine frequentiel 
• hypothese : reponse periodique 
• outil : HBM, methode de balance harmonique 



_ t Nh 

series de Fourier : U(t) = U 0 +^ t/„. c cos(/7CO/) + 0'„. i sin(/70V) 

n=l 



I 

procedure de Galerkin: f(U r ) — Z r U r + X— F r 

multiplicateurs de Lagrange 
representant les forces de contact 



• nouvelle expression du probleme complet 
• en respectant les contraintes NL definissant 
• les forces de contact exprimees dans le domaine temporel 
• c’est un probleme d’optimisation avec contraintes 



• les ingredients de la strategic 

• aspect non-lineaire : le processus est iteratif, solveur de type Newton 
• aspect probleme contraint en dynamique : Dynamic Lagrangian [Nacivet, 2003] 

• souci temporel / frequentiel : alternating Frequency Time [Cameron et Griffin, 1989] 



S. Nacivet, these ECL, 2003 L. Salles, these ECL, 2010 
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Amortissement / challenge recherche 



Comment I’amortissement structural evolue-t-il avec I’usure ? 



Pression 
Aerodynamique 



Nceuds de 
contact 




Force 

centrifuge 



Noeud 

d’excitation 

harmonique 




Evolution de I’amplitude 
vibratoire au sommet de 
I’aube au cours du temps, 
avec usure o.e- 



I’amplitude augmente 



I’amplitude diminue 




o.i 

0.S95 



1 .005 



Normalized Frequency 



Normalized fretting cycle 




L. Salles, these ECL, 2010 
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Amortissement / challenge recherche 



Comment departager I’effet amortissant d’un dispositif de I’effet du fluide ? 




D. Laxalde, these ECL, 2007 
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Prise en compte du couplage aeroelastique 

Line demarche 

• La presence du fluide modifie les caracteristiques modales estimees in vacuo; 
reciproquement le mouvement d’une aube influe sur les mouvements du fluide 

• un calcul couple est done a priori necessaire, en particulier en presence de non-linearites 

• une hypothese consiste a realiser un couplage faible, permettant des calculs separes pour 

le fluide et la structure, I’effet du fluide est alors double : 

1/ effet d’amortissement aero ; 2/ effet de forgage 



STRUCTURAL UNSTEADY AERODYNAMICS STFA AERQDVNAMiCS 




Schematic of forced response prediction system [Hilbert et al. 1997] 
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Prise en compte du couplage aeroelastique 

Estimation de I’amortissement aero 



• selection des modes d’interet (etude des coincidences, importance des dephasages inter-aubes) 
• calcul aero avec mouvement modal impose des aubes 
• obtention des forces de pression parietale instationnaires 
• plusieurs voies sont alors possibles : 

Voie A 

• calcul du travail des forces de pression dans le mouvement de I’aube pour un cycle 

* estimation du ratio Travail cyclique / Energie cinetique moyenne 
• analogie amortissement modal materiau / aero via le concept de decrement log par exemple 

Voie B 

* incorporation des forces aero dans le calcul force via des matrices d’amortissement 

proportionnelles aux deplacements de I’aube 




Damping (%) 


Euler Nun Linear 


NS linearized 


Structural 


0.17 


aero elastic 


0.26 


0.31 


Total 


(1.41 


0,47 



Table 1 : damping estimation 



M. Jocker, these KTH, 2002 / E. P. Petrov, GT2009-59634 / exemple : Seinturier etal, GT2004-53372 
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Prise en compte du couplage aeroelastique 

Calcul de la fonction de forgage 



• les pressions parietales instationnaires sont des forces imposees sur le modele structure 
- on s’interesse a la premiere harmonique et eventuellement aux suivantes 

- les charges sont def Mantes 

• resolution numerique du probleme : cf cours F. Thouverez 




Figure 9: Unsteady NS total pressure (1* L harmonic 
modulus) on pressure (PS) and suction sides (SS) 




Figure 11: Comparison of predictedttest results 
(blade dynamic levels vs. rotating speed) 

NB : modele avec mistuning (recale) 



exemple : Seinturier etal, GT2004-53372 






Plan 



A. Pre-design aero 

B. Analyse du comportement sous chargement statique 

C. Estimation du risque vibratoire 

I. Calcul des modes propres 

II. Determination des frequences d’excitation 

III. Recherche des coincidences 

IV. Calcul des contraintes dynamiques 

V. Trace du diagramme de Haigh 

D. Tests sur prototypes 




C.V. Trace du diagramme de Haigh 



Principe 



• notion de limite d’endurance 



Conlraimt 
k njpi ure 



S* 

iimlle 

il'pn^urnnr^ 




Tflfl 1 QP 1 Q 2 |fl 3 



10 * tO 5 IQ 1 10 * 10 ^ 

Nrambr* de CVCteS A fUptiJrtf 



FEgurrs Q « Lns trn-iji dcmnin^s d'tindunnce 




MO 



N (cycles) 



1 

8 



j £50 



Traction 



,p= 9 '^ 



itssesen TA6V 



• effet de la contrainte moyenne sur I’amplitude 
de contrainte alternee maxi admissible 



Contrainte 

fillKrnEB 




T | i d r ; 3 i l pa GtiOdm&n 3 : fiaritLiOle de GErber 

2 : drake de Sk'debery 4 : construction VDI 

Figure- T3 - piflgr-frmme* de Hfllgh 



Techniques de I’lngenieur, BM 4 175 et 5 052 / LAMEFIP Bordeaux 





Plan 




A. Pre-design aero 

B. Analyse du comportement sous chargement statique 

C. Estimation du risque vibratoire 

D. Tests sur prototypes 

I. Extensometrie 

II. Vibrometrie par laser Doppler 

III. Mesures par tip-timing 





D.lll. Tests / Extensometrie 



Principe 




Encyclopedia of vibrations. Elsevier, 2002 / Test devices Inc. 







D.lll. Tests / Extensometrie 



Caracterisation vibratoire sous chargement aero 
des redresseurs du banc CREATE (1/3) 





RM1 

(64) 



RM2 RM3 

1 I 



RD1 



RD3 






D.MI. Tests / LDV 



Principe 






Wikipedia : laser vibrometer/ www.polytec.com 






D.lll. Tests / Tip-timing 



Principe 




• Les sondes sont montees sur les parties fixes. 

• On mesure les temps d’arrivee des aubes devant chaque sonde. 

• Les vibrations provoquent des ecarts entre temps mesures et calcules a partir des tops arbre. 

• On en deduit les frequences et les amplitudes des modes d’aubes. 





Conclusions 



* des mecanismes multi-physiques et multi- 
echelles 

• des couplages de plus en plus marques : 
technologies monobloc, allegement, charge / 
aube accrue 



in 




Rot:K Uva speed 

Figure- 20: predictionftest comparison - gage 2 



• un fort besoin de mieux comprendre les 
phenomenes pour de meilleures pratiques de 
conception 




FUTURE project (http://www.future-project.eu), analyses performed by 
Volvo Aero (http://www.volvoaero.com), Siemens Industrial Turbomachinery AB 

(http://www.sit-ab.se), 

MTU Aero Engines GmbH (http://www.mtu.de), DLR - Deutsches Zentrum fur Luft- 
und Raumfahrte.V (http://www.dlr.de/), 

PCA Engineers Limited (http://www.pcaeng.co.uk), ONERA (http://www.onera.fr/), 
T urbomeca (http://www.turbomeca.com/), 

Techspace Aero (http://www.techspace-aero.be/ ), Snecma 
(http://www.snecma.com/) 



Seinturier etal, GT2004-53372 / FUTURE resultats intermediaires 






Perspectives 




• couplage solidien arbres BP-HP / roues aubagees : nature des modeles, NL [M. Gruin] 

• couplage fluide-structure fort + desaccordage + amortissement [G. Benefice] 

• systemes dissipatifs actifs [B. Zhou] 

• travaux theoriques et experimentaux sur le pompage, le flottement ... 

• et aussi : les excitations transitoires, les SFV ... 
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